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A key element in the architecture of a quantum information processing network is a reliable phys- 
ical interface between fields and qubits. We study a process of entanglement transfer engineering, 
where two remote qubits respectively interact with entangled two-mode continuous variable (CV) 
field. We quantify the entanglement induced in the qubit state at the expenses of the loss of entan- 
glement in the CV system. We discuss the range of mixed entangled states which can be obtained 
with this set-up. Furthermore, we suggest a protocol to determine the residual entangling power of 
the light fields, inferring, thus, the entanglement left in the field modes which, after the interaction, 
are no longer in a Gaussian state. Two different set-ups are proposed: a cavity-QED system and 
an interface between superconducting qubits and field modes. We address in details the practical 
difficulties inherent in these two proposals, showing that the latter is promising under many aspects. 

PACS numbers: 03.67.-a, 03.67.Mn, 42.50.Pq, 85.25.Dq, 85.40.-e 
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I. INTRODUCTION 

Distributed networks for quantum communication and 
quantum computation have recently received large at- 
tention being a promising architecture for quantum in- 
formation processing. Efficient sources of entangled con- 
tinuous variable (CV) states of light are readily available 
and field modes can act as reliable information carriers. 
On the other hand, it is more handy to manipulate the 
information stored in qubits embodied in atomic or solid- 
state systems. In many of the protocols for quantum in- 
formation processing designed for distributed networks, 
entanglement plays a crucial role in establishing an ex- 
ploitable quantum channel between two distant nodes . 
This motivates the attempt to study and implement a 
reliable physical interface between fields and qubits. An 
interface allows for the transfer of information from the 
carrier to the qubit subsystem. In this context, it is in- 
teresting to consider the ability of a quantum-correlated 
two-mode field to induce entanglement in several differ- 
ent qubit subsystems. 

On the other hand, under a more fundamental point of 
view, an analysis of the entangling power of a correlated 
two-mode state can be a useful tool to indirectly quantify 
the entanglement left between the modes after the inter- 
action with a given pair of qubits. For a non-Gaussian 
field state, in particular, there is a lack of objective crite- 
ria to determine whether or not entanglement is present 
between two modes. We show that the capability of the 
light fields to induce entanglement in two initially sepa- 
rable qubits provides a test for the inseparability of the 
fields. 

The implementation of such a physical interface, thus, 
opens a way to investigate the exchange of entanglement 
between systems defined in heterodimensional Hilbert 
spaces 0: S and to the entanglement transfer pro- 



cesses, where entanglement increases in one system at 
the expense of the loss of the other one, through their 
interaction. 

In this work, we study a two-qubit system interacting 
with a quantum-correlated two-mode squeezed state of 
light through bi-local resonant interactions. The system 
we analyze allows to engineer entanglement transfer, gen- 
erating qubit states with a controllable amount of entan- 
glement. We show that it is actually possible to explore 
a large part of the space of the entangled mixed states 
(EMS) of two qubits, including the class of boundary 
states having the maximum amount of entanglement for 
a given value of the purity of the state @ . The generation 
of arbitrary and controllable entangled mixed states is 
relevant to investigate about the role that entanglement 
and purity have for the tasks of Quantum Information 
Processing (QIP). It has been proved, for example, that 
while entanglement is a fundamental requirement to effi- 
ciently perform the Shor's factorization, purity is not |^. 
Moreover, in some cases, efficient QIP can be performed 
even relaxing the requirements, in terms of quantum cor- 
relations and purity, of a quantum channel. For example, 
a neat threshold exists for the allowed mixedness of an 
entangled mixed state to be usefully used in a teleporta- 
tion protocol 0- 

The interest in the generation of EMS and in quanti- 
fying the entangling power of a given two-mode entan- 
gler leads us to look for practical implementations of the 
model we propose here. We consider two different set-ups 
to implement our model, namely a cavity-quantum elec- 
trodynamics (cavity-QED) system and Josephson charge 
qubits interfaced to field modes. We describe the two 
systems in some details and compare their performances. 
In the first proposal we consider two optical cavities, each 
interacting with one of the modes of a two-mode squeezed 
state of light respectively. After the squeezed state feeds 
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the cavities, two two-level atoms cross their respective 
cavities. However, this set-up faces important practical 
difficulties. The cavity-photon lifetime turns out to be 
comparable to the operation times of our set-up, a fea- 
ture that is typical of the weak coupling regime of the 
atom-field interaction. On the other hand, the advan- 
tages of a much stronger coupling j^, can be exploited 
in the second scenario we propose, where the qubits and 
the cavity are implemented by superconducting devices 
integrated on the same chip. In this case, flying qubits 
are no more required so that the control over the in- 
teraction times is easier. Moreover, in this scheme the 
qubit parameters can be independently modulated via 
external electric and magnetic fields The interac- 

tion can be controlled by tuning the qubit on/off res- 
onance with a cavity mode |3j 13 ■ A promising design 
of integrated superconducting qubit and cavity has been 
recently proposed [U and the first experiments have al- 
ready demonstrated a quality factor Q > 10'* for the 
cavity which is enough to implement the protocols 
discussed in refs. 0, . 

The paper is organized as follows. In Section^we de- 
scribe the general scheme proposed here, introduce the 
space of the entangled mixed states of two qubits and 
show that our model allows for an extensive exploration 
of such a space. In Section ITTll we infer the correlations 
left between two field modes after the interaction. We 
propose a scheme that, extending the interaction model 
to more than a single qubit pair, allows to get some addi- 
tional information on the residual entanglement capabili- 
ties of the field, bypassing the lack of a necessary and suf- 
ficient criterion for the entanglement of a non-Gaussian 
state. In Section IIVI we describe the proposed imple- 
mentations of our system. For a cavity-QED system, we 
discuss the effect o f sp ectator atoms and of the cooper- 
ativity parameter This kind of difficulties are no 

more present when Cooper pair boxes [loj | integrated in 
high-quality transmission lines are used. 



II. THE MODEL OF ENTANGLEMENT 
GENERATION 

We introduce the system we consider, based on the 
interaction of a pair of remote qubits with local environ- 
ments that, for definiteness, we model as single bosonic 
modes, respectively. While the physical system that em- 
bodies the qubits will not be specified until necessary, this 
schematization perfectly matches the situation in which 
the two qubits are coupled to the field modes of two cav- 
ities. This assumption does not limit the generality of 
our approach allowing, on the other hand, to cover many 
different physical systems that are promising for the pur- 
poses of QIP 01 . We will consider explicitly this case, 
from now on. We are interested in a situation in which 
the field modes of the two cavities exhibit non-classical 
correlations and local interactions of each qubit with the 
respective cavity mode are arranged. A way in which this 



can be achieved is sketched in Fig. ^ 
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FIG. 1: Two remote cavities are driven by a quantum cor- 
related two-mode field. We explicitly consider the case of a 
two-mode squeezed driving field coupled to each cavity via a 
leaky mirror and preparing a quantum-correlated two-mode 
state of the cavity fields. After the preparation of the cav- 
ity field, local interactions with two generic (flying or static) 
qubits are arranged. 

In details: two remote but identical cavities are ini- 
tially prepared in their vacuum states. One mirror of 
each cavity is assumed to be perfect while the other 
has non-zero transmittivity. The leaky mirror is cou- 
pled to one mode of an external two-mode squeezed state 
with coupling strength k. The squeezed state of the 
two external modes a and b is represented by 15')^^ = 
(coshr)"^ J2'^,oi^^^^'^T l")a \^)b' where {\n)}^j^ are the 
Fock bases for the field modes, and r is the squeezing 
parameter. If the coherence time of the driving field is 
shorter than k~^, we can treat the cavity mirror as a 
beam splitter continuously fed by squeezed fields [l5l |. 
The state of the cavity modes A and B evolves due to 
the coupling and is described by the reduced density ma- 
trix ^2] 



PAB = Tr,b (BAa{O)BBb{O)pabAB{O)B\,{90Bbi^) 



OO min(n,m) 

= E E I" - ~k\(g> \n - l)j,{m - l\ . 

n,m— k,l—0 

(1) 

Here, BAa{0) — exp 9{Ad^ — A'^a) is the beam split- 
ter operator 0| with its refiection coefficient sin 61 (ba- 
sically determined by k I15|), d (d^) is the annihilation 
(creation) operator for the external field mode a (analo- 
gously for field mode b) while A (A^) and B (B^) are the 
bosonic operators of the cavity fields. We have defined 
XriniG^l^, with Xnm = (tauh r)"+™/(cosh r)^ and 
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When the refiected part of the external field is minimized, 
a two-mode squeezed state quickly builds up inside the 
cavities. Note that a two-mode Gaussian field is built 
into the cavity, regardless of the coupling strength. This 
model is exactly the same as the initial field interacting 
with a zero-temperature reservoir [Tsj and it is known 
that entanglement can always be found during the evo- 
lution of an initial two-mode squeezed field in the zero- 
temperature vacuum |T^ . 

After the cavity field is prepared, the interaction with 
a pair of qubits with logical states |0) ■ and |1) - {i — 1, 2) 
begins. Here, the interaction model is assumed to be of 
the resonant Jaynes-Cummings type T^. The interac- 
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tion Hamiltonian for the cavity A is 

HAi = hn (i|i)i(o| + it|o)^(i 



(2) 



where is the atom- field couphng constant. An analo- 
gous Hamiltonian describes the interaction between the 
second qubit and cavity B. Before entering the detail of 
protocols, we point out that the Hamiltonian we dis- 
cuss can be also implemented by superconducting de- 
vices mum. 

The interaction between cavity modes and qubits gives 
rise to an unitary evolution Uai (t) ®Ub2 (t) of the whole 
system, where tJAi{t) = exp{—iHAit/h). The effective 
evolution of the two qubits, on the other hand, is non- 
unitary and is described by the reduced density matrix 
Pi2{t) obtained by tracing out the cavity fields as pi2(t) = 

Ttab (uAi{t) (g> t>B2(i)pi2(0) (S> PAsUiiit) (g> J/laW) ■ 

Here, we consider the initial state pi2(0) = |00)j2(00| 
and define the rescaled time of this first interaction 
n = nt. In the basis {|11) , |10) , |01) , |00)}i2, pMn) 
takes the form 



P12(ti) 



/ A -D\ 
B 
C 
\-D OOF 



(3) 



with F = 1- A- B - C, 

OC n 

A = J2Y1 -4fcrsin^(ri\/;r^)sin2(TiV^r^), 



n=0 k,l=0 

Qo n 



B = ^ ^ A^" sm^{TiVn- k) cos^(TiVn - I), (4) 

11=0 k,i=a 

oo n 

C = J2J2 ^rcos2(nV^r^)sin2(ri\^r^), 

n=0 k,l=0 

and the off-diagonal component given by 
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D = J2J2 -^k?^^ smirWn- k + l) 

X cos(TiVn — k) sin(ri-\/n — ^ + 1) cos{TiVri^^) . 

(5) 

In order to quantify the quantum correlations between 
the qubits after the interaction with the cavity modes, we 
use the negativity of the partially transposed density op- 
erator (NPT), which is a necessary and sufficient condi- 
tion for entanglement of any bipartite qubit system |2C|. 
The entanglement measure based on NPT is defined as 
^NPT = ~2A~, where A~ is the negative eigenvalue of 
the partially transposed density matrix (here with 
respect to qubit 2) ^21]. In our case, A~ does not depend 
on the populations of states |ll)i2 ; |00)i2- Explicitly 



In Fig. [51 we plot £n pt versus ti for a reflectivity of 
the cavity sin^ 6* = 0.1 and three different values of the 
squeezing parameter r. The entanglement is peaked at 
Ti = (2g -I- l)7r/2 {q integer). For small values of r, 
the two-mode squeezed state can be approximated by 
|00)^^ + ''|H)a6 ^^^^ ^^"^ qubit-field interaction re- 
sults in |00)^2 ~ ^ |ll)i2 '''1 ~ (2? + l)'''/2- However, 
as squeezing is increased, the Rabi oscillations become 
more complicated due to the importance of terms rela- 
tive to higher photon numbers in eqs. 10)) and ((SJ). This 
reduces the qubit entanglement and makes £npt a non- 
monotonous function of the squeezing. This is explicitly 
shown in Fig. [2 for r — 1.2 (full squares). In this case, 
the entanglement function is always below the curve rel- 
ative to the smaller r — 0.86 (stars) that allows for the 
maximum £npt at ri w 37r/2. 

Another useful quantity to characterize the state of 
the qubits is the degree of purity of pi2 . We have stud- 
ied the behavior of the linearized entropy Sl{pi2) = 
(4/3) [l — Tr (P12)] , which gives a measure of the de- 
gree of mixedness for the state described by pi2 
Sl{pi2) ranges from (pure states) to 1 (maximally 
mixed states). The linearized entropy can be studied 
as a function of ri and for different values of r. We 
have combined the temporal behavior of entanglement 
and linear entropy onto an entanglement-purity plane as 
shown in Fig. O In this plot, each curve is relative to 
a specific value of the squeezing parameter, while each 
point along a curve gives the entanglement and purity 
of the corresponding state at a specific interaction time 
(so that Ti represents a curvilinear abscissa, in this plot). 
The solid line is the upper bound to the region occupied 
by physically achievable mixed states for a given degree 
of entanglement Enpt > 0. They are known as maxi- 
mally entangled mixed states (MEMS) . A parameter- 
ization of MEMS is critically dependent on the chosen 
measures of entanglement and purity. If the negativity 



•'NPT 




X- = ]^\b + C - ^ AD^ + [B - Cf 



(6) 



FIG. 2: Enpt versus the rescaled interaction time r (in units 
of tt). The squeezing parameter of the field is r = 0.26 (full 
line), r — 0.86 (stars), r — 1.2 (full squares). These values 
offer a trade-off between the visibility of each curve and the 
insight given in the dynamics of the entanglement. We take 
sin^6l = 0.1. 
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FIG. 3: Navigation in the plane of the entangled mixed states 
of two qubits. We used sin^ 6 = 0.1. The measures for the en- 
tanglement and purity are negativity of partial transposition 
Enpt and linearized entropy Sl- The curvilinear abscissa is 
ri G [0, 37r/2]. We show the cases with r = 0.46 (rhombuses), 
r — 0.86 (stars) and r — 2 (filled squares). The solid line is 
the MEMS boundary. 



of partial transposition is taken, there are two classes of 
one-parameter states belonging to MEMS and giving the 
same boundary. The first class corresponds to the family 
of Werner states pw ~ p\4'^){4>^\ + (1 ~ p)^2x2, with 
10+) = {1/V2) (|00) -I- 111)) and < p < 1. The other 
class can be resented by 
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(7) 



i±4±2Z; 



Thus, the Werner states represent a frontier for the 
achievable amount of entanglement that a given mixed 
state can have. This feature is unique for the choice of 
NPT as the measure of entanglement [5| . 

Several interesting points can be addressed closely 
analysing Fig. O First of all, the interaction model de- 
scribed here is able to produce entangled, nearly pure, 
state of two qubits (see, for example, the last point on 
the curve relative to r = 0.46) that can be used to test 
protocols for QIP Q. As the value of r increases, the 
states of the qubits get closer to the frontier curve of 
MEMS, even if the corresponding states are weakly en- 
tangled for longer periods of time. As was pointed out, 
this is due to the contribution of highly excited photon 
number states, in the Schmidt decomposition of a two- 
mode squeezed state \S), for larger values of r Our 
theoretical model is flexible enough to generate arbitrary 
entangled mixed states of two qubits up to the boundary- 
class of MEMS. Recently, theoretical and experimental 
efforts have been performed in the generation of these 
states using different mechanisms |22l |. In our scheme. 



the local interactions with the two-mode entangler in- 
duces time-dependent quantum correlations between the 
qubits. The mixedness of the two-level systems state is 
due to their entanglement with the cavity modes. We will 
see that, in fact, the entanglement between the qubits is 
set at the expenses of the correlations between the field 
modes. 



III. ENTANGLEMENT DYNAMICS FOR 
CAVITY FIELDS 

In the previous section, we studied the amount of en- 
tanglement generated in a bipartite two-level system by 
the interaction with a two-mode squeezed state. Now, 
we turn our attention to the entanglement left between 
the cavity fields. After the interaction, the cavity fields 
state is given by 

PAB{t) = Tri2 {UAB12{t)pAB{0) ® Pl2 (0)?7^bi2 W) , (») 

where lJABi2{t) = UAi{i) ® UB2{t)- It is not difhcult to 
see that the cavity field is no more Gaussian after the 
interaction. 

Differently from the qubit case, how to quantify the 
quantum correlations of a CV state is only partially 
known. For the case of a Gaussian state, NPT can be 
used as a separability criterion as well as a measure 
of entanglement [23. 12^ . In this case, the NPT con- 
dition for separability is equivalent to the violation of 
the uncertainty principle by a partially transposed den- 
sity operator (23j. To study the uncertainty principle, 
it is useful to consider the vector of the field quadra- 
tures I = (9i,pi,g2,P2)^, where qi = (A + A'') /V2, 

Pi 

fined in terms of B and i?^. The field quadratures satisfy 
the commutation relations [Cq,^/?] = i^ap, where flap 
are the elements of the 4x4 matrix ft = 0^J, with 
1' 



j4 — j4M / \/2 and 52 , P2 are analogously de- 



J = 



Some of the statistical properties of a 



two-mode CV state can be inferred from the 4x4 covari- 
ance matrixV, defined by Vap — (A^q,A^^-|-A^^A^q)/2, 

with A^^ = ia — (ia^ ^^'^ ^^'^ expectation values eval- 
uated over the state of the light fields. In terms of V, 
the uncertainty relation for the field quadratures takes 
the form V + ift > 0. Using the block representation 
A 



V = 



B 



with A, B and C 2 x 2 matrices, the 



uncertainty relation can be restated as |22 



; (det A)(detB) 
(det A + detB) 



> 0. 



detC)^-Tr(ACBC'^) 



(9) 



with A — AJ, B — BJ and C — CJ. Each term of 
A is invariant under local linear canonical transforma- 
tions. Partial transposition is equivalent to a mirror re- 
flection (that performs the transformation pi —pi) in 
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the phase-space. This changes the sign of det C, leaving 
the other terms unaffected Thus the uncertainty 

^NS for the partially transposed density matrix is ob- 
tained replacing C by |det C| in Eq. ©. In our study, the 
function A^s depends on the squeezing parameter r and 
the interaction time ti. For a Gaussian state, Ans < 
is a necessary and sufficient condition for entanglement. 

For non-Gaussian states the situation is more com- 
plicated. The violation of the Heisenberg uncertainty 
relation by is only a sufficient condition for non- 
separability of the state. Nevertheless, as far as the au- 
thors are concerned, the uncertainty criterion is one of 
the most successful conditions to test entanglement of 
even a non-Gaussian state. We challenge this condition 
in the following. 

In Fig. 01 Ans for the cavity fields after the interaction 
with the qubits is plotted against ri for two significant 
values of the squeezing parameter. As done before, we as- 
sume the qubits initially prepared in their ground states 
|00)j^2- Interestingly, Ans is negative for most of the in- 
teraction time, when r = 0.46 (solid line). This means 
that the cavity field modes are in an entangled state even 
after the qubit-field interaction. The dynamics of A^vs is 
opposite to the qubit entanglement as it is apparent com- 
paring Figs.[51and01 This behavior is consistent with the 
idea that the entanglement in the field is transferred to 
the qubits. On the other hand, for r = 0.86 the entangle- 
ment induced between the qubits does not really match 
to the behavior of Ans (dashed line). 



that is the maximum achievable entanglement £npt,m 
for a given interaction time ti. We consider r — 0.86 
and all the qubits initially in their ground state. Obvi- 
ously, if the first pair of qubits does not interact with 
the cavities (ti — 0), the entanglement settled in the sec- 
ond pair is the maximum achievable. Comparing Figs. [21 
and [51 we note that the second pair may be entangled 
more as the first pair is entangled less and vice-versa. 
Notice that at t ~ 1.5, i.e. when the entanglement of the 
first pair is maximal (see Fig. j^J, we find that £npt,34 
is nonzero, showing that even in this case the entangle- 
ment capability of the cavity fields is not exhausted by 
the first interaction. Thus the field modes in the cavi- 
ties are still quantum-mechanically correlated and able 
to entangle the two qubits by the local interactions, de- 
spite the fact that Ans is positive. In this example, the 
entanglement capability is a more powerful test for the 
entanglement of the cavity fields. 

maxxj£NPT,34 




^NS 
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FIG. 4: This plot shows the dynamics of the uncertainty 
relation with partially transposed density matrix. Negative 
values show the violation of uncertainty or, equivalently, the 
non-separability of the cavity field modes A and B depending 
on the atom-field interaction time. We set (sin ef =0.1 with 
r = 0.46 (solid line) and r = 0.86 (dotted line). 

More informations about the entanglement of the field 
modes can be obtained testing their entangling power, 
that is their ability to set entanglement in an additional 
pair of initially separable qubits. In this case, if the cav- 
ity fields are not entangled there is no way to transfer 
correlations to two additional qubits. A nonzero entan- 
gling power is, thus, a condition for the fields states to 
be inseparable. 

We let a second pair of qubits (indexed by 3 and 4) 
interact with the respective cavities for a time T2. In 
general, the degree of entanglement £npt,3.a between the 
qubits depends also upon ri, the interaction time with 
the first pair. In Fig. we plot max^-^ 5ArpT,34('''i, T2), 



FIG. 5: This plot shows the entanglement between qubits 3 
and 4, maximized with respect to the interaction time T2, as 
a function of ri (units of vr). Because of the limit in compu- 
tational power, the maximum of £npt,34, is calculated for the 
cases of the first pair interaction times ri = nn/4 (n=0,l,..,8). 



IV. TWO PHYSICAL IMPLEMENTATIONS 

The generation of two-qubit quantum states up to the 
MEMS boundary and the interest in inferring the entan- 
gling capabilities of a non-Gaussian CV state motivate 
the search for an implementation of the model studied 
so far. In this Section we describe two different propos- 
als for a set-up: a cavity-QED scheme and an interface 
between superconducting charge qubits and cavity field 
modes. This latter, in particular, offers some intriguing 
perspectives in terms of coherence times and control of 
the qubits. 



A. Cavity-QED system 

The first physical set-up we analyze is sketched in 
Fig. O As outlined in Section m two remote one-sided 
optical cavities are initially prepared in the vacuum state. 
A two-mode squeezed state is coupled to each cavity via 
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the leaky mirrors. In what follows, we assume a broad- 
hand external source. We call Awe^t the band-width of 
the dyriving source. To consider the external field as an 
infinite band-width drive and use the (simple) analytical 
expressions valid for a broad-band squeezed state, the 
condition n ^ Awe^t has to be fulfilled. However, it is 
typically lS.ujext ^ ^ [T^ . This means that to compare 
any theoretical prediction to the results of a real experi- 
ment, a more involved finite band- width approach must 
be used . This goes beyond the tasks of this work and 
we will assume the broad-band condition for our theoreti- 
cal investigation. We assume the cavity fields build up as 
a two-mode squeezed state before the interactions with 
the qubits start. 

The qubits are here embodied by two fiying two-level 
atoms of their ground and excited states \g)^ , |e)j (i — 
1, 2) that pass through the respective cavities. 
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FIG. 6: Two remote optical cavities are fed by the two-mode 
squeezed state generated by a nondegenerate-parametric- 
amplifier (NOPA) . The external fields interact with two cavity 
modes via direct coupling to the leaky mirror of each cav- 
ity (the other being perfect). A low-density flux of two-level 
atoms, initially prepared in their ground states \g), passes 
through each cavity and interacts with the relevant cavity 
mode. 

The behavior of a one-sided cavity, with respect to an 
external driving field, is influenced by the presence of an 
atomic medium. A two-level atom in a resonator can be 
seen as an intra-cavity lossy medium with a loss parame- 
ter proportional to the atomic cooperativity C = fl"^ / kT, 
where F is the spontaneous emission rate from |e). For 
large C, the intra-cavity losses are so large that, eventu- 
ally, no squeezed state builds up in the resonators [l3l| . 
Realistic parameters, within the state of the art, are 
(k, fl, F)/27r = (80, 40, 4) MHz. While the interaction of 
a cavity mode with the respective flying atom arises natu- 
rally in the form of a Jaynes-Cummings Hamiltonian, an 
important issue to take into account is the strength of this 
interaction. The optimal condition would be, obviously, 
the strong coupling regime 17 3> k, F, where the coherent 
evolution of the atom-cavity system is faster than the 
decoherence mechanisms (dissipation, dephasing) due to 
the cavity decay and to the atomic spontaneous emis- 
sion [2^. In general strong coupling with flying atoms 
is a hard task to achieve, in optical systems 26]. Fur- 
thermore, with the above values for k, suitable for the 
external field to be injected, the field coherence times 



are within the range of lOnsec. On the other hand, 
for cavity field waist of ~ 10 fj,m and an atomic velocity 
~ 300 m/ sec (that is an optimistic value) we get transit 
times in the range of 10 nsec, comparable to k~^. The 
dissipative effect due to the cavity decay is, thus, not 
negligible and experimental efforts are required to give 
a practical implementation of our scheme in this set-up. 
This has some implications for our method to infer the 
entangling capabilities of the fields by detecting the cor- 
relations between the atoms of a second pair. Intuitively, 
the time elapsing between the passages of two consecutive 
atoms has to be shorter than the life-time of a photon in 
an optical cavity. Otherwise, the decay of the cavity fields 
will destroy the quantum correlations between the field 
modes. The effects of the cavity decay can be described 
by introducing a dissipative term Hioss = —ifiK'^^a^a 
{a = A, B) in the system Hamiltonian. This term gives 
an exponential decay (with rate k) of the probability to 
flnd the cavity mode a in a Fock state with p photons. 
We call t the time elapsing between the passage of the 
flrst and the second pair of qubits. It turns out that, for 
i equal to an atomic transit time (so that the first pair 
of atoms has surely left the region of interaction with 
the cavities), £npt,3A is 50% less than what we get in 
the ideal conditions. A way to minimize the elapsing 
time is to have the simultaneous presence of at least two 
qubits in the same cavity. For the sake of definiteness, 
we suppose the intensity of each cavity field to have a 
Gaussian radial profile centred at the cavity axis. We 
assume that, while the first atom is interacting with the 
cavity field, the second is at the border of the region of 
interaction and is weakly coupled to the field. We thus 
have a spectator atom inside the cavity |l3l | . Usually, the 
spectator atoms give rise to additional loss mechanisms 
that become important once the density of the spectators 
is such that their collective coupling to the cavity mode 
is comparable to Q.. A way to reduce these losses would 
be, then, to control a very low-density atomic beam. 

Having fixed the interaction time (that can be finely 
controlled), a possible source of errors, in our model, is 
represented by the mismatched injection of the atoms 
in the cavity. The distribution of the atomic velocities is 
thermal and it is, in general, hard to arrange the simulta- 
neous entrance of two atoms in two remote cavities. Ob- 
viously, a certain control is possible by means of atomic 
cooling techniques. However, a quantitative analysis of 
the effect of a mismatched triggering of the qubit-cavity 
mode interactions shows that the qubit-qubit entangle- 
ment is flattened to zero for times short compared to the 
interaction time ti. The long-time behavior of £npt.12, 
then, follows the pattern expected for a perfectly injected 
pair of qubits. In general, for a delay St between the en- 
trance of 1 and 2 much smaller than il^^, this effect 
on the generation of a pair of entangled atoms is neg- 
ligible. Obviously, a complete description of the effect 
of a mismatched atomic injection is given averaging the 
entanglement over St. Such a detailed analysis is not 
necessary whenever we are able to keep the delay times 
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within the condition St ^ fJ"^. However, the control of 
dt is per se a hard task. Some of the problems faced by 
this optical cavity-QED implementation can be solved in 
a scenario in which microwave cavities, interacting with 
long-lived Rydberg atoms, are considered. In this case, 
however, there is still the need for a fine control of the 
transit times and the simultaneous entrance of the atoms 
inside the cavities. 



B. Josephson qubits in a superconducting 
transmission line 



We now consider a second class of implementations 
that combines some of the features of a microwave cavity- 
QED with the characteristics of a superconducting nano- 
circuit. The first advantage of a solid state device com- 
bined with quantum optics is the possibility of achieving 
a strong qubit-field mode Jaynes-Cummings type inter- 
action Q . This is basically due to the fact that the cav- 
ity field, in this case, is coupled to the charge of the 
qubit rather than its dipole. Various solid state imple- 
mentations of the cavity have been proposed, ranging 
from large Josephson junctions |MI3to superconducting 
films with large kinetic inductance to microstrip res- 
onators Several operational schemes for quantum 
protocols have been suggested, in this contexts 0, 0,113 ■ 
Couplingbetween a superconducting qubit in the charge 
regime 0, and a classical Josephson junction has 
been implemented by the Saclay group [23. Decoher- 
ence of a qubit coupled with a quantum oscillator in the 
solid state has been studied in ref. This analysis has 
revealed the existence of optimal operating conditions 
where decoherence is essentially due to spontaneous non- 
radiative decay of the qubit and leakage from the com- 
putational space of the cavity modes. This latter source 
of decoherence is minimized if the cavity is im plem ented 
by an integrated high-Q microstrip resonator so, in 
what follows, we will focus on this particular implemen- 
tation. 

In the proposal of ref. 'ill, ^ Superconducting- 
Quantum-Interference-Device (SQUID) qubit 10] is 
placed between the planes of the resonator and the 
whole device is fabricated using nanolithographic tech- 
niques. The geometry of the resonator is such that a 
single mode of frequency 10 GHz can be accommo- 
dated. Microstrip transmission-lines having a Q fac- 
tor of about 10'' (corresponding to photon life-time of 
0.1/isec) have been already built and the effect of cou- 
pling with a qubit implemented by a SQUID has been 
demonstrated |1^]. Realistic estimates of the qubit life- 
time are in the range of some fisec and the whole sit- 
uation is reminding of a microwave cavity-QED set-up 
where the internal losses due to the cavity decay are very 
low. The source of quantum radiation can be built up 
using nondegenerate Josephson parametric amplifiers (in 
a distributed configuration, to limit the effects of gain- 
increasing noise). Microwave squeezed radiation (in a 



range of frequency of ~ 10 GHz) has been demonstrated 
and the source impedance-matched to superconducting 
lines used to transport the signals |30l |. 

A sketch of the set-up we consider is given in Fig. [7| 
(a). Let us first concentrate on the interaction of a 
SQUID qubit with a single cavity mode. We operate 
at the SQUID degeneracy point where the qubit is en- 
coded in equally-weighted superpositions of states hav- 
ing zero and one excess Cooper-pair on the SQUID is- 
land, namely |±) = (1/^2) (|0) ± |2e)) (2e being the 
charge of a Cooper pair). The free Hamiltonian of the 
SQUID is given by Hsquid = ^Ej{(j)ext)^z- Here, is 
the z-Pauli operator and Ej((j)ext) is the Josephson en- 
ergy which is tuned via an external magnetic flux (pe^t 
piercing the SQUID loop. This changes the energy sep- 
aration between |-|-) and |— ) and can be used to switch 
on/off resonance the interaction with the field mode. The 
degeneracy point is set by biasing the SQUID with a dc 
electric field via the ground plate of the resonator (Fig.[7| 
(a)). To quantify the strength of the qubit-field interac- 
tion, we model the cavity mode as a LC oscillator coupled 
to the SQUID as in Fig.[7|(b). The coupling is realized 
via the capacitor Gc- The interaction Hamiltonian can be 
cast in the form of a Jaynes-Cummings interaction with a 
Rabi frequency f2 that, for proper choices of the circuital 
parameters, is as large as 0.5 GHz 0, 0|- The strong 
coupling regime is thus possible in this set-up and inter- 
action times of about ~ lOnsec ^ k^^,T^^ are enough 
to generate an entangled state of two qubits near to the 
MEMS boundary. 



(a) 



_ ^ capacitive gap - 



I qubit 1 




microstrip resonatorA 



FIG. 7: (a): Coupling between the field mode of a mi- 
crostrip resonator and two SQUID qubits. We show just one 
resonator-qubits subsystem. This scheme has to be doubled 
considering a resonator B interacting with qubits 2 and 4. 
Qubit 1 and 2 are used to generate a MEMS, while the inter- 
action of the cavity modes with qubit 3 and 4 is useful to infer 
the state of the field modes, (b): Quantum circuit showing 
the LC-oscillator model (parameters Lq, Co) for the cavity 
mode. The coupling with the qubit is capacitive (coupling 
capacitance Cc). Cjg is the capacitance of each Josephson 
junction. 

Two SQUID qubits (size ~ fim) can be easily accom- 
modated in the cavity (L ~ 1 cm) far enough to achieve 
negligible cross-talking (in principle due to direct capac- 
itive and inductive coupling). Lithographic techniques 
allow to control (within few percent) the geometric char- 
acteristics and the resulting parameters of the device. 
The two qubits can be manipulated both simultaneously 



8 



with a uniform magnetic field or independently with two 
separate coils. Due to charged impurities in the vicin- 
ity of the device, separate calibration at the degeneracy 
points is required for each qubit. This may be achieved 
with several adjustments of the design of ref. for in- 
stance by splitting the ground plate and by attaching a 
gate to each part. 

The qubits have fixed position inside the cavities. 
This is an advantage with respect to common microwave 
cavity-QED systems, where flying atoms are required. 
The interaction times are regulated tuning the interac- 
tion between the SQUID and the cavity mode on and off 
resonance via (/)ext- This avoids the problem of different 
qubit injections in the cavities even if the two SQUIDs 
have to be set on resonance with the respective cavity 
mode at the same time. The results we derived in Sec- 
tion about a mismatched interaction-triggering are 
fully applicable here. On the other hand, having pre- 
cisely fixed the number of qubits inside the cavities, we 
do not have to deal with spectator qubits. 

We now briefly address the effect of decoherence in 
this set-up. Working at the charge degeneracy point, we 
minimize the effect of charge-coupled sources of noise . 
At this point, decoherence due to low-frequency modes 
vanishes at first order, a key property which allows to 
achieve coherence times of several hundreds of nanosec- 
onds in the Saclay qubit p9l |. Furthermore, in our set- 
up, a selection rule prevents direct transitions between 
the states of the dressed doublet 0. As a result of this 
analysis, the performance-limiting process turns out to 
be the spontaneous decay from each state of the dressed 
doublet to the ground state. This results from compa- 
rable contributions of non-radiative decay of the qubit 
due to 1// noise and losses due to the resonator. Life- 
times of ^ 300 ns were estimated for a system where 
the resonator is implemented by a Niobium Josephson 
junction. As we said, the losses of the resonator are fur- 
ther minimized in the proposal of ref. 11], where 1/f 
noise is the ultimate limiting factor. Noise sources are 
switching charged impurities and their effect depends on 
statistical properties of the environment which are be- 
yond the power spectrum [sij. The actual dephasing 
rate depends also on details of the protocol and may 
show device-dependent features |32|]. A detailed analy- 
sis with simulations of noisy gates for this device will be 
presented elsewhere "ss]. It is worth stressing here that, 
due to the qubit-resonator interaction, the energy levels 
of our qubit are much less sensitive to charge fluctuations 
than an isolated qubit at the optimal working point. This 
leads to a conservative estimate of coherence times in the 
range of 100 ns [s^, which surely allows for navigation 
and generation of MEMS. 

At the degeneracy point, computational states of the 
qubit cannot be distinguished by measuring their charge. 
To detect the state of a SQUID qubit, we have to slowly 
shift the working point far from the degeneracy point, 
sweeping a dc-bias, to adiabatically transform the states 
of the qubit as \+) |2e) , |-) |0). Then, charge 



measurements can be performed |28l |. 

As far as the interaction with more than one pair of 
qubits is concerned, we make use of the further degree of 
freedom represented by the tunability of the energy spac- 
ing between the levels of each qubit. We can proceed as 
follows. We suppose SQUID 1 and 3 in resonator A (as in 
Fig. □(a)) while SQUID 2 and 4 interact with B. First, 
just 1 and 2 are resonant with the corresponding cavity 
mode while 3 and 4 are in a dispersive regime obtained ei- 
ther dc-biasing them or using an external magnetic flux. 
Once the interaction time ti has passed, we set the inter- 
action with this pair in a dispersive regime while 3 and 
4 are set on resonance with the respective cavity mode. 
The timing of these operation can be controlled electron- 
ically and the operating time can be as large as 50 nsec. 
In this condition, an entangled state of 3 and 4 is estab- 
lished still within the coherence time of the system. 

In summary this second set-up offers some advantages 
with respect to a cavity-QED implementation. The most 
important points are related to the longer coherence 
times of the dynamical evolution while this is the main 
limitation of an optical set-up. An important problem 
for this solid state system is 1/f noise. In our opinion 
further improvements can be achieved by development of 
the design exploiting the possibilities of nanolithography. 



V. CONCLUSIONS 

When a quantum-correlated CV state, such as a two- 
mode squeezed state of light, interacts with a bipartite 
two-dimensional system via bi-local interaction, effective 
entanglement transfer is possible. Theoretically, the state 
of the qubits can be engineered in terms of entanglement 
and purity. The model provides a tunable source of en- 
tangled mixed states that can be useful to investigate 
the interplay between entanglement and purity in QIP 
and for purposes of quantum communication and com- 
putation. The entangler, i.e. the CV system, does not 
exhaust its entangling capabilities with a single interac- 
tion but is able to entangle other pairs of qubits. This 
property can be exploited as a test, based on the entan- 
gling power of the field, for the quantum correlations in 
a non-Gaussian state of light. 

We have proposed two different set-up in which this 
scheme can be implemented. The first is a cavity-QED 
system in which the qubits are embodied in two-level 
atoms crossing two optical cavities. The second pro- 
posal exploits the recent ideas about solid-state sys- 
tems/quantum optics interfaces and uses superconduct- 
ing qubits integrated in microstrip resonators. This sec- 
ond scenario, in particular, offers the advantages of a 
strong coupling regime of interaction (that is hard to get 
with optical cavities) without the difficulties connected 
with the management of flying qubits. 
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